Using All-Sky Surveys to Find Planetary Transits 



Joshua Pepper, Andrew Gould and D. L. DePoy 

Department of Astronomy, Ohio State University, Columbus, Ohio 43210 

pepper@astronomy . ohio-state . edu , gould@astronomy . ohio-state . edu , 
depoyOastronomy . ohio-state . edu 

O 

(N 



O 

Q 

o 

in 



ABSTRACT 

Transits of bright stars offer a unique opportunity to study detailed prop- 
erties of extrasolar planets that cannot be determined through radial-velocity 
observations. We propose a new technique to find such systems using all-sky 



>- ■ small-aperture transit surveys. We derive a general formula for the number of 

^3 ■ stars that can be probed for such systems as a function of the characteristics of 

the star, the planet, and the survey. We use this formula to derive the optimal 

00 ' 

telescope design for finding transits of bright stars: a 5 cm "telescope" with a 
g ■ 4k x 4k camera. 

Subject headings: techniques: photometric - surveys - planetary systems 
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<3 ■ 1. INTRODUCTION 

> '■ 

^ . In the past three years, great strides have been made in the detection of extrasolar 

planets (XSPs). To date, nearly all of the roughly 100 known XSPs have been discovered 
using the radial velocity (RV) technique. However, RV detections, in and of themselves, yield 
only a few planetary parameters, namely the period P, the eccentricity e, and M sin(i), where 
M is the mass of the planet and i is the inclination of its orbit. By contrast, if a planet transits 
its host star, much more information is available. First, of course, the Msin(z) degeneracy 
can be broken. Second, the ratio of the radii of the planet and host star can be measured. 
Therefore, provided that the star can be classified well enough to determine its mass and 
radius, then the planet's radius and hence its density can be determined. Third, and perhaps 
most important, if the transits can be observed with sufficient signal-to-noise ratio (S/N), 
then one can probe otherwise unobservable details of a planet, such as its oblateness (Hui 
& Seager 2002), atmospheric conditions (Charbonneau et al. 2002), and perhaps satellites 
and rings. Regardless of how a planet is initially discovered, once it is determined to transit 
its host star, this wealth of information can in principle be extracted by intensive follow-up 
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observation of these transits. This fact has been amply demonstrated by the discovery and 
analysis of the transiting planet HD209458b (Charbonneau et al. 2000; Cody & Sasselov 
2002). 

At the moment, all ongoing and proposed transit surveys are carried out in relatively 
narrow pencil beams. They make up for their small angular area with relatively deep ex- 
posures. These surveys fall into two basic classes: field stars (Howell et al. 2000; Brown & 
Charbonneau 1999; Mallen-Ornelas et al. 2001; Udalski et al. 2002) 1 , and clusters (Street et 
al. 2000; Burke et al. 2002) 2 . These surveys are potentially capable of establishing the fre- 
quency of planets in various environments, but they are unlikely to find the kinds of transits 
of bright stars that would be most useful for intensive follow-up analysis. Although some 
of the surveys of field stars are considered "wide field", their total survey areas are small 
compared to An str. One project that has the potential to cover a very large area is WASP 
(Street et al. 2002), which plans to employ five cameras, each with a 9.°5 x 9.°5 field of view. 

An alternative method is to conduct an all-sky survey. Instead of continuous observation 
of all targets (which is impossible from a practical standpoint for an all-sky survey), this 
approach would necessarily involve revisiting each target in the sky at regular, semi-regular, 
or random intervals throughout the course of the project. This kind of observing strategy will 
not yield a continuous light curve on any star, as the current transit surveys do. Rather, this 
plan will generate an only sporadically sampled light curve. However, the long time baseline 
for the survey will eventually generate just as many individual observations of a single star. 
Transit-like dips in the data stream will not be visually obvious, but by repeatedly phase- 
folding the full light curve back on itself over a range of periods, one can detect the dips 
from the transits. (See § 3.) 

This approach is especially relevant given the fact that there are several all-sky surveys 
already being planned for objectives other than transit detections. It should be possible, 
for instance, to utilize the photometric data stream of upcoming astrometric missions for 
transit detection. Space-based projects such as GALA 3 and DIVA 4 would take hundreds 
of observations of millions of stars over mission lengths of years with the aim of obtaining 
precise astrometry. These data could equally well be analyzed for planetary transits. 



1 http: / / www.psi.edu / ~ esquerdo /asp / asp. html 
http://www.hao.ucar.edu/public/research/stare/stare.html 
http : / /bulge . princeton.edu / ~ ogle / 

2 http://star-www.st-and.ac.uk/^yt2/WEB_GROUP/top.html 

3 http: / / astro.estec.esa.nl/GAIA/gaia.html 

4 http://www. ari.uni-heidelberg.de/diva/diva. html 
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There are also several existing or proposed ground-based all-sky surveys, including the 
Large-aperture Synoptic Survey Telescope (LSST) 5 , the Panoramic Optical Imager (POI), 
and the All-Sky Automated Survey (ASAS) 6 (Pojmanski 2000). These surveys will be imag- 
ing the entire sky every few days, a qualitatively similar cadence to those of GAIA and DIVA. 
While LSST and POI will likely be saturated by stars of V < 12, which are of greatest in- 
terest for transit follow-ups, ASAS is of particular interest in the present context because of 
its very small aperture. 

In this paper, we examine the process of analyzing photometric data streams from all- 
sky surveys to find transits. We calculate the sensitivity to XSP detection, the distance out 
to which these detections will be possible, and the number of false-positive detections due to 
random noise. We derive a general expression for the number of stars that can be probed by 
this technique as a function of the total sensitivity of the survey. We apply our analysis to 
the problem of telescope design and conclude that very small, 5 cm telescopes are optimal 
for finding transits of bright stars. 



2. SCALING RELATIONS 

In order to accurately determine the number of transits a given survey would be expected 
to detect, we must carefully define the set of stellar systems that can be probed for transits 
by the survey. This will evidently depend on the number density n, the luminosity L, and 
the radius R, of the stars being probed; as well as the semi-major axis a and radius r of 
their planets. A quick (and naive) formulation would then state that the total number of 
systems that can be probed for transits N p is 

4 R 

N p = -Tin— [d max (L,R,a,r)f (1) 

o (X 

for a homogeneous population of such stars, where d max is the distance out to which a transit 
can be detected. 

However, the quantity d max is ill-defined. For fixed L, R, a, r and distance d, the transit 
detection scales as (1 — x 2 ) 1 ! 2 where x is the transit impact parameter normalized to R; that 
is, < x < 1. Hence, the signal-to-noise ratio scales as (1 — x 2 ) 1 ^. Since detections normally 
require a minimum S/N, d max must also be a function of the impact parameter x. Thus, for 



5 lit t p : / / w w w . lssto.org/lssto/index. htm 
6 http://www. astrouw.edu.pl/~gp/asas/asas. html 
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a photon- noise limited survey, d max oc (1 — x 2 ) 1//4 , and therefore 

-y 2 - = t n U ~ \ d raax (L, R, d, T, x)f (2) 

ax 6 a 

where d max (L, R, a, r, x) = d max (L, R, a, r, 0)(1— rr 2 ) 1 / 4 . Here d max (L, R, a, r, 0) is the distance 
out to which a transit can be detected for an edge-on (i = 90°) orbit. We must then integrate 
over all values of the impact parameter x from to 1, and so 

P = J ~dV = -^^maA^R^^^f (3) 

where, 



n= [\l-x 2 f 4 dx = y^^§-^ 0.719. 
'Jo 2 T(9/4) 



(4) 



We now determine the dependence of N p , the total number of systems probed, on the re- 
maining parameters L, R, a, and r. To do so, we analyze the detection requirement, 

N t (^j >A X 2 min , (5) 

where N t is the number of observations of the transit over the length of the survey, 5 is the 
fractional change in the star's brightness during the transit, a is the fractional error of an 
individual flux measurement, and Ax 2 ^ is the minimum acceptable difference in x 2 between 
a fit that assumes a constant flux and one that takes account of a transit. As we discuss in 
§ 3, Ax 2 lin must be set sufficiently high to avoid spurious detections due to random noise. 



• To determine the dependence of N p on L, we note that in equation (3), the only factor 
that depends on L is d^ nax . For a particular star, the flux / = Lj (4nd 2 ), thus d oc L 1 / 2 , 
and so N p oc L 3/2 . 

• For the dependence of N p on R, we note that / oc d~ 2 , and that a oc f' 1 ^ 2 , so d oc a. 

1 /2 

From equation (5), we see that a oc 5N t . Since Nt is the total number of observations 
of the star during transits over the length of the survey, N t = N y }S (2R)/(27ra). Also, 
5 = (rrr) 2 /(ttR) 2 . Thus, d max oc R~ 3 / 2 . Combining these relations with the explicit 
factor of R in equation (3) itself, we arrive at N p oc Rr 1 ! 2 . 

• For the dependence of N p on a, we see that d oc a, and a oc 5N^ 2 . Using N t oc a -1 , 
we have c^ aa . oc a~ 3//2 , and so N p oc a~ 5 / 2 . 

• Finally, for the dependence of N p on r, the only factor that depends on r is 5, from 
5 = {jir) 2 /(TiR) 2 . Thus a oc r 2 , and so N p oc r 6 . 
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Consolidating the dependence of N p on the various parameters, we finally arrive at, 

where do = d max (L Q , Rq, a , r , 0). That is, d is the distance out to which a planet-star 
system with i = 90° and the fiducial parameters L , Rq, ao, r can just barely be detected at 
the S/N threshold. 

We now seek to simplify equation (6) by integrating over the local stellar population at 
fixed absolute magnitude, My, and so we replace the three independent variables (n, R, L) by 
the single variable My. We consider two regimes My > 6 and My < 6 (with one overlapping 
bin at My = 6, which we will use later to check for consistency). 

We first treat the fainter regime. Here, the main sequence is relatively narrow. Hence, 
R may be regarded as a function of L (and so, therefore, of My), while n is simply the 
number density of stars in a given magnitude bin. Hence, the "integration" amounts to a 
simple multiplication of factors. 

We adapt the number density of stars n from the empirically determined local stellar 
luminosity function (LF): for the range (9 < My < 18) we use the LF reported in Zheng et 
al. (2001), and for the range (6 < My < 8) we use the LF of Bessell & Stringfellow (1993). 

To estimate stellar radii, we combine the linear color-magnitude relation, My = 3.37(V— 
J)+2.89 from Reid (1991), a color/surface-brightness relation \og(R/R Q ) = 0.69 + 0.2226(V- 
K) — 0.2M V , based on the data of van Belle (1999), and VIK color-color relations for dwarfs 
from Bessell & Brett (1988). We can therefore calculate the relative number of systems with 
a fixed a and r as a function of My, which we designate 



F(M- 
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n(M } 



v 



n 



L(M V 



3/2 



R{My) 
Rq 



-7/2 



(7) 



where n ,L , and Rq are the normalizations chosen below. 



For the upper main sequence, My < 6, we evaluate F(My) directly using the Hipparcos 
catalog (ESA 1997). For example, the LF for My = 4 would be computed by summing 
^ i [(4/3)7rDf] _1 over all stars within the Hipparcos completeness limit, V < 7.3, having 
3.5 < My < 4.5, and lying within 50 pc. The distance A is the minimum of 50 pc and 
the distance at which the star would have V = 7.3. Actually, we are not directly interested 
in the LF at My = 4, but rather in the integral of V'^Rr 1 ! 2 over the subpopulations that 
make up the My = 4 bin of the LF. Hence 

3.5<M Vii <4.5 u u 
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The stellar radii are determined from Hipparcos/Tycho (B T , Vt) photometry and the color /surf ace- 
brightness relation of Gould & Morgan (2003), \ogR/R Q = 0.597 + 0.536(5 T - V T ) -M Vt /5, 
ultimately derived from van Belle (1999). 

To normalize F, we adopt the values associated with My = 5 stars: L = 0.86 L Q , 
Ro = 0.97 Rq, n = 0.0025 pc" 3 . 

The resulting function (Fig. 1) shows that the majority of stars that are probed will 
be F and G type (2 < M v < 6.5). To check how this distribution depends on the size 
of the volume sampled, we recalculate the distribution function for the case in which the 
observed stars cover a much larger volume - one which would be better described by a thin 
disk, rather than a spherically uniform distribution. For this case we find that the scalings 
shown in equation (7) are replaced by F(My) oc nhLR~ 2 , where h is the scale height of 
each population. The distribution function is still dominated by F and G stars, but there 
are more K and early to mid M stars (7 < My < 12). Of course, it is common knowledge 
that magnitude-limited (F oc nL 3 ^ 2 ) samples of main-sequence stars will be dominated by 
F and G stars. The interesting feature of Figure 1 is that this result does not qualitatively 
change despite the addition of the factor i?~ 7//2 in equation (7), which very strongly favors 
later-type stars. 



3. RANDOM NOISE 

Equation (6) describes what kinds of XSP systems can be detected by a certain survey, 
given a photometric detection limit. The threshold Ax^ in is determined by taking account 
of the fact that the data stream from an XSP search must be analyzed for any combination 
of the parameters R, a, and r within reasonable ranges. Such analysis will, however, yield a 
number of false-positive detections due to random noise. The threshold value of Ax^ in must 
be chosen to yield a manageable number of candidate systems for follow-up observations. 

To determine Ax^ in , we generate 1000 independent simulated streams of photometric 
observations of a single system with a host star of one solar mass and a circular planetary 
orbit. We attempt to simulate a schedule that would be characteristic of an all-sky survey 
that re-images a given star approximately every few days with varying intervals between 
observations. For each data stream, we generate 1000 observations at irregular intervals over 
1800 simulated days, and then phase the observations for a range of periods from 3.0 to 3.1 
days and a range of transit lengths. The number of the observations and the duration of 
the simulated survey are arbitrarily chosen as plausible characteristics for the type of all-sky 
survey we envision. 



-7- 



The duration of the transit depends on the orientation of the system with respect to our 
line of sight, and we analyze the data for eight equally spaced, progressively more inclined 
orientations. For each orientation, we test for 16 equally spaced phases. We test the resulting 
data sets for transit-like dips in the light curve, which we define simply as intervals during 
which the local mean light curve dips significantly below the global average. This (8 x 16) 
grid structure is chosen empirically: we find that the number of false positives increases 
linearly with grid density below this density and then flattens above it. 

The result shows that for this kind of system, in order to restrict follow-up analysis to 
the 0.1% of the full sample most likely to yield a true transit detection, the value for Ax^ in 
should be set to ~ 36.6. We take the highest value of A% 2 from each of the 1000 sets, and 
of those highest values, we then sort the 1000 sets from highest value of A% 2 to lowest. In 
Fig. 2 (inset), we plot the highest value of A% 2 for each of the 1000 data streams. 

To check the robustness of this number, we run additional simulations of 25 data sets, 
each with different observation schedules. First, we generate a set of observations that 
are randomly distributed throughout a 1800 day mission, with the requirement that no 
two observations be less than 10 minutes apart. Then we regenerate the data set with 
observations that are evenly spaced throughout the project. We also conduct the analysis 
on data sets with evenly spaced pairs of observations, with the observations in each pair 
separated by 14.4 minutes; and then again with the pair-spacing at 43.2 minutes. 

These different configurations test for different types of observing schedules. For in- 
stance, a ground-based all-sky survey would most likely observe a star at essentially random 
times throughout a project. On the other hand, a space-based mission with a slowly rotating 
telescope (similar to the Hipparcos satellite) would observe a star in pairs of observations 
separated by minutes. For various reasons there could be a certain regularity imposed on 
either the space-based or ground-based observations. 

We find that the value for A^ 2 ^ does not depend strongly on the observing schedule. 
We determine this by taking the highest value of A% 2 for each of the 25 independent data 
streams, sorting the highest values of A% 2 from each set in the manner described above, and 
then plotting the results for each type of observing schedule (Fig. 2). We find that there is 
only a ~ 10% variation in Ax 2 ^ for different schedules. 

The value of Ax 2 ^ depends on the number of different parameter configurations that are 
tested. When we rerun our analysis using a period range of 3.0 to 3.5 days, we increase the 
size of our parameter space by a factor of 5. We expect that the value of Ax 2 ^ depends on the 
size of the parameter space, according to Ax 2 min (^ new) - AxLmi^oid) = 2 ln(^ new /^ oM ). 
This prediction agrees with the simulations, which show that an increase in parameter space 
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by a factor of 5 leads to an increase in Ax^ in by 21n(5) = 3.2. 

We then predict ^Xmm ^ the size of the parameter space is expanded to search for 
transits with periods between 2 and 10 days (the range of periods most probable for the sort 
of fiducial values of the other parameters we have chosen). This would increase \1> by a factor 
of 80 compared to the period range of 3.0 to 3.1 days. Therefore, Ax^ in for the expanded 
range of periods is higher than 36.6 by 21n(80) = 8.8. So the value expected for Ax^ in for 
such an analysis is about 45. 



4. NUMBER OF SYSTEMS PROBED 

The weak dependence of Axm in on the observing schedule implies that the sensitivity of 
a project to planets essentially depends only on the total number of photons detected from 
each star, and not on the details of how they are collected. For a given star, this number is 
obviously proportional to the flux. We therefore characterize the sensitivity of the observing 
setup (telescope + detectors + duration + weather + etc.) by 7, the total number of photons 
that are detected from a fiducial V = 10 mag star during the entire project. (Here we adopt 
V = 10 as a reference, although the stars of interest lie in the range 8 < V < 10.) We can 
utilize the various relations used to derive equation (6) to relate 7 to V max , the maximum 
apparent magnitude at which an equatorial transit can be detected, 

V„» = -2.51o g (^k^!) +10 . (9) 
Then, considering equations (5), (6), and (7), we obtain 



(10) 



where we have adopted 70 = 1.25 x 10 7 , a = 10 R Q , r = 0.10 R Q , and where we have made 
our evaluation at M v = 5 (i.e. R = 0.97 R Q , V max = 10, d = 100 pc, and n = 0.0025 pc -3 ). 
Note that 70 = 1.25 x 10 7 corresponds to approximately 625 20-second exposures with a 5 
cm telescope and a broadened (V + R) type filter for one V = 10 mag fiducial star. 

As mentioned above, a noteworthy feature of equation (10) is that N p depends on the 
characteristics of the survey primarily through the parameter 7. Moreover, since Ax^ in 
depends only logarithmically on the size of the parameter space being explored, it plays a 
minimal role in survey design compared to the other variables in equation (10). 



We envision two scenarios to which these results will be applicable. In one, a stream of 
photometric measurements from a space-based astrometric mission, such as GAIA, could be 
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searched for transits. In the other, a ground-based survey could use one or more dedicated 
telescopes to search all bright stars in the solar neighborhood. 

5. IMPLICATIONS FOR TELESCOPE DESIGN 

We now apply the general analysis of § 2 and § 4 to the problem of optimizing telescope 
design for quickly locating a "large" number of bright (V < 10) transiting systems. Since 
only one such systems is now known, we define "large" as (9(10). From equation (10) and 
the 0.75% frequency of hot jupiters measured from RV surveys, there are roughly 5 such 
systems to be discovered over the whole sky per magnitude bin at My = 5 for V max = 10. 
Hence, from Figure 1, of order 25 are to be discovered from all spectral types. It would, of 
course, be possible to discover even more by going fainter, but setting this relatively bright 
limit is advisable for three reasons. First, as we argued in the introduction, the brightest 
transits are the most interesting scientifically, and most of the transits detected in any survey 
will be close to the magnitude limit. Second, as we discuss below, a high dynamic range, 
AV = — V mm , can only be achieved at considerable cost to the observing efficiency 
Hence, if high efficiency is to be maintained, setting V max fainter means eliminating the 
brightest (most interesting) systems. Third, at V max = 10, we are already reaching distances 
of 100 pc for G stars. Hence the number of transits observed in fainter surveys will not 
continue to grow as djj as in equation (6). 

In previous sections, we ignored the loss of sensitivity to systems that are brighter than 
V min , which is set by saturation of the detector (or more precisely, by the flux at which 
detector non-linearities can no longer be accurately calibrated). This fraction is 10~ 0SAV , or 
6% for AV = 2, which we therefore adopt as a sensible goal. That is, we wish to optimize 
the telescope design for, 

8 = V min < V < V max = 10. (11) 

(In any event, essentially all stars V < 8 have already been surveyed for XSPs using RV, and 
the problem of determining which among the planet-bearers have transits is trivial compared 
to the problem of conducting an all-sky photometric variability survey.) 

Optimization means maximizing the photon collection rate, 7/T, where T is the duration 
of the experiment and 7 is, again, the total number of photons collected from a fiducial 
V = 10 mag star. Explicitly, 

7 = KSTD 2 ^^, (12) 

where A8 is the angular size of the detector, D is the diameter of the primary-optic, S 
is the fraction of the time actually spent exposing, and K is a constant that depends on 
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the telescope, filter, and detector throughput. For our calculations, we assume K = K = 
40 e~ cm -2 s _1 , which is appropriate for a broad (V + R) filter and the fiducial V = 10 mag 
star. The design problems are brought into sharper relief by noting that A6 = C/DJ 7 , where 
C is the linear size of the detector and T is the focal ratio, or //#, of the optics. Equation 
(12) then becomes 

KSC 2 T 

7 = "4^- (13) 

That is, almost regardless of other characteristics of the system, the camera should be 
made as fast as possible. We will adopt T = 1.8, below which it is substantially more difficult 
to fabricate optics. A more remarkable feature of equation (13) is that all explicit dependence 
on the size of the primary optic has vanished: a 1 cm telescope and an 8 m telescope would 
appear equally good! Actually, as we now show, there is a hidden dependence of £ on D, 
which favors small telescopes. 



5.1. Considerations for Aperture Size 

The global efficiency £ can be broken down into two factors, £ = £o£s, where £ is the 
fraction of time available for observing (i.e., during which the sky is dark, the weather is 
good, etc.), and £ s is the fraction of this available observing time that the shutter is actually 
open. The first factor is not affected by telescope design and so will be ignored for the 
moment. The second factor should be maximized. The smaller the telescope aperture is, the 
longer the exposures can be before a V m i n = 8 mag star saturates. Since the readout time 
is fixed, a smaller fraction of time is lost to read-out. We adopt as benchmarks a detector 
with pixel size of Ax p = 9 /zm, well depths of 10 5 e~, and a telescope diameter of 5 cm. 

To make explicit calculations, #psf, the full width at half maximum of the point spread 
function (PSF) must be specified. For the fast optics {T < few) we consider here, the 
diffraction limit is always much smaller than a pixel, regardless of aperture: 0^/9, p ~ 
1.22 JFA/Axp ~ 0.16 (for our fiducial choices). At the small apertures we will consider, the 
diffraction limit is larger than the seeing, so it is possible to make the PSF much smaller 
than a pixel, #psf <C 6 p . This would have the advantage of reducing sky noise and is a 
useful approach when it is possible to always center the telescope at the same field position 
as is the case for "point and stare" experiments. However, for an all-sky survey, which 
cycles through many fields, such precision repeat pointing is extremely difficult. Without it, 
precision photometry is impossible unless the sub-pixel response of the CCD is mapped out in 
detail. We therefore adopt a Nyquist-sampled PSF, for which the sky noise is approximately 
that falling on 4n ~ 13 pixels. 
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Our overall consideration for telescope design must take into account three factors. 
First, we with to maximize observing efficiency Second, we wish to achieve the highest 
possible signal-to-noise ratio. Third, we must avoid any distortion problems with the optics. 
There are four effects through which aperture size can impact these factors. Two of these 
effects, observing efficiency and scintillation noise, will drive us to larger telescopes, while the 
other two effects, sky noise and focal plane distortion, will drive us to smaller telescopes. As 
we show below, for the observing parameters we have specified an aperture of 5 cm ensures 
a manageable (and unique) balance between the various effects. 



5.1.1. Exposure Time vs. Readout Time 

Assuming Nyquist sampling, at most half the light from a point source falls within one 
pixel. We can directly calculate the ratio of time lost to readout T rca d to the time spent 
exposing T cxp , 

T r ead , ( D \ n -QA(v„„-M) ( W \ K T r 



T 

± exp 



1 in-0.4(V ma x-10) K Tread 1 n0.4(Ay-2) / -i a\ 

\5cmJ \WJ K 30s ' 1 ' 



where W is the well depth of the detector pixels, and W = 10 5 e~ is a fiducial well depth. 
Note that the factor lo -0 - 4 ^™ -10 ), which arises from the need to avoid saturation of the 
brightest stars (where V m i n = V max — AV), has been broken up into two terms to permit easy 
comparisons of equation (14) with equations (16) and (18) below. 

In order to maximize the efficiency £s, the fraction of observing time devoted to readout 
should be minimized, and therefore, according to equation (14), so should the aperture size. 
The telescope will operate reasonably efficiently so long as T rcad < T cxp . 



5.1.2. Scintillation 

Another concern that arises for small apertures is the effect of atmospheric scintillation, 
which is characterized by (Young 1967; Warner 1988), 

^^J^Y m x^^(-^)(^Y'\ (15) 
1 \cm J \ h J \ sec / 

where Sq = 0.09, X is the airmass, h is the altitude of the observatory, and ho = 8 km is 
the scale height of the atmosphere. Using values of D = 5 cm, X = 1.5, h = 2 km, and 
T cxp = 30 s, we find AI/I = 0.0057. 
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Hence, for exposure times set by a saturation threshold, T cxp oc D 2 (see eq. [14]), 

Scintillation Noise / D 
Source Noise \5cm 

Therefore, despite the common perception that scintillation is a greater problem for smaller 
telescopes, for the fixed photon counts per exposure that are of interest in the present 
context, scintillation noise increases with increasing aperture. However, this dependence is 
fairly weak. For D = 5 cm, the scintillation noise is just slightly smaller than the photon 
noise at V = 10. 



I/O 



10 



-0.2(V-10) 



(16) 



5.1.3. Sky Noise 



All the calculations in §§ 2, 3, and 4 have assumed that sky noise is negligible, i.e. 
V <C V s k y where V^k y is the light from the sky falling on 13 pixels (for Nyquist sampling). 
Assuming a somewhat conservative mean sky brightness of V = 20.0 mag/arcsec 2 , one finds 



V sky = 10.6 + 5 log 



D T {Ax. 



5 cm 1.8 \9 /im 



-i 



;i7) 



To determine how much of a problem sky noise will be, we can compare it to the amount of 
photon noise, 

Sky Noise 3 / D 
Source Noise 4 V 5 cm 



-i 



10 0.2(V-10) 



(18) 



Thus, sky noise is less than photon noise for a 5 cm telescope, but would become a serious 
problem for a substantially smaller aperture. 



5.1.4- Focal Plane Distortion 

Focal plane distortions toward the edge of the detector become difficult (and expensive) 
to correct when the field of view is too large. For example, for a Ak x Ak detector with 
Ax p = 9 /xm pixels, 

A ^ 23 °(^)"(n)"- (19 » 

Again, this is manageable for D = 5 cm but could would potentially be a problem for smaller 
apertures. 



Although this field of view is not so large as to create focal plane distortions, it is 
important to note that for a field of view this large, the telescope must be placed on an 
equatorial mount. For the alternative, an alt-az mount, the rotation of the sky will cause 
stars at one edge of the field to move faster across the detector than stars at the opposite 
edge. For similar reasons, the telescope must track rather than using drift scan. 



5.2. Optimal Telescope Design 

The relationships described in § 5.1.1 through § 5.1.4 can be used to optimize the 
aperture D for any survey parameterized by a given V max . One is driven to smaller apertures 
by the goals of minimizing scintillation noise and the fraction of time spent on readout, and 
to larger apertures by the goals of minimizing sky noise and field distortion. Given available 
C = 3.6 cm 4096 x 4096 detectors and reasonably fast T = 1.8 optics, a D = 5 cm telescope 
is optimal for an all-sky survey of V = 10 stars. Among all existing transit programs of 
which we are aware, the WASP telescope (D = 6 cm lens, T = 2.8 focal ratio, 2k x 2k, 
C = 3 cm detector, Street et al. 2002) comes closest to meeting these design specifications. 

We use equation (13) to determine the required duration of the experiment using our 
optimally-designed telescope, adopting 7 = 70, K = 40e~cm~ 2 s~ 1 , £ = 20%, £3 — 50%, 
C = 3.6 cm, and T = 1.8. The total time required to conduct the survey assuming only 
source photon noise is T = 4 months. Taking into account sky noise and scintillation increases 
that time by a factor of 2.5 to 10 months. 

For surveys that intend to search for transits of stars fainter than the range we consider 
in this paper (8 = V min < V < V max = 10), the relations in § 5.1.1 through § 5.1.4 operate 
somewhat differently. For the magnitude range we are considering, equations (14) and (16) 
provide an upper limit on the aperture size, while equations (18) and (19) place a lower 
limit, in which both limits converge to D — 5cm. At fainter magnitudes, the aperture size 
is restricted according to equations (14) and (18): 

1 q0.2(V-10) < -P < ]n0.2(V max -10) 

5 cm 

This equation shows how the aperture size limits combine into a single scaling relation. What 
of the other two limiting factors? The issue of focal plane distortions is unimportant at fainter 
magnitudes, since A6 < 23°, which is outside of the regime where such distortions are a 
factor. Scintillation effects are also not a factor at fainter magnitudes, since the scintillation- 
noise restriction requires that (_D/5cm) < io°- 6 ( l/_10 ) ) which is a looser restriction than that 
of equation (20). 



(20) 
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Therefore, for a survey of transits at a magnitude range 8 = Kiin <V< V max = 10, the 
aperture size limits converge to D — 5 cm, while at fainter magnitudes, the aperture size is 
given by equation (20). 

5.3. Practical Implementation 

In the previous section, we estimated the duration of observations required to achieve 
the minimum S/N to detect transits by hot Jupiters assuming certain fiducial parameters 
of a ground-based telescope, but calculated within the framework of a literal "all-sky" (47r) 
survey of randomly-timed observations that is more characteristic of satellite missions. The 
resulting estimate is useful for judging the viability of a given observing setup, but it glosses 
over a key issue in the detection of transits, namely the problem of folding the data. As 
discussed in §3, the number of folds (and hence the size of the search space that must be 
probed) scales directly as the duration of the experiment. This larger search space increases 
both the minimum A% 2 for a robust detection and the amount of computing power needed to 
sift through the search space. The first effect is logarithmic in the size of the search space, so a 
10-fold increase changes AxLm by only ~ In 10 ~ 2.3, which is well under 10%. However, the 
second effect is linear in the search space and so could easily overwhelm available computing 
resources if not carefully monitored. That is, there are important drivers for keeping the 
duration of observations of any given field to a minimum and hence for exploring the question 
of whether it is better to break up the "all-sky" survey into several smaller components, each 
of which could be completed in a shorter time. Indeed, the OGLE experiment (Udalski et 
al. 2002), the only transit experiment to successfully detect a transit (Konacki et al. 2003), 
was motivated by these considerations to concentrate its observations over a month duration 
so as to limit the number of foldings. 

To estimate the true duration of the project required the achieve the minimum S/N, one 
must take account of two factors: First, as a practical matter, ground-based surveys from a 
single location can cover only an angular area Q < 4ir. Second, during a year of continuous 
observations, any given patch of sky is observed only for about 6 months. These factors 
change the previous estimates from §5.2. The first factor implies that the observing time 
required to reach minimum S/N is actually lower than the time calculated from equation 
(13) by a factor of Q/4n ~ 0.5, since the project is observing about half the angular area 
on the sky. The impact of the second factor is more complex. In a given night, only half 
the accessible angular area Q can be observed. Therefore, each night the available area to 
observe is lower by an additional factor of 2, which means the rate of observations of a given 
point on the sky will double again. There are two scenarios at this point. The first scenario 
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is that after accounting for the first factor, the time required for sufficient S/N is less than 
a year. In this case, the two factors combine to lower the estimate from §5.2 by a factor of 
about 4. In the second scenario, the time required for sufficient S/N is more than a year. In 
this case, after observing for 6 months the minimum S/N is not reached, and therefore the 
project will have to pick up again six more months later when the target is again visible, 
and so the second factor does not apply. As we see below, however, for the parameters we 
are considering, we are well within the regime of the first scenario. 

Applying the first factor to the result in §5.2 reduces the calculated time from 10 months 
to 5 months. Since this latter duration is indeed smaller than 1 year, the second factor implies 
that the requisite S/N will actually be reached in 2.5 months. This is only slightly longer 
than the duration of the OGLE observations. 

However, real experiments inevitably have larger errors than expected. (We mention one 
source of additional errors below.) If the errors prove sufficiently larger that the experiment 
requires more than 1 year (after application of the first factor) to achieve the minimum S/N, 
then one would not gain the advantage of the second factor. In this case, it would be better 
to break the sky up into strips by declination, and observe each strip for a year, so as to 
increase the S/N obtained in each strip during a single year, and so to permit the application 
of the second factor. 

Another real-world consideration is that the photometric errors will not be Gaussian. 
For example, the Sloan Digital Sky Survey (SDSS) photometry errors, while Gaussian in 
their core, deteriorate to an exponential profile in the wings beginning at about 3 a (Ivezic, 
Z. et al. 2003). Such deterioration is likely to set in earlier for the small-aperture, wide- 
field, low-budget cameras that we envisage here. The non-Gaussian form of the errors has 
no practical impact on our analysis: there are so many data points that the central limit 
theorem guarantees that their combined behavior in each phase bin will be Gaussian (as we 
have implicitly assumed in §3. However, the non-Gaussian tails will tend to increase the a 
of the distribution relative to what would be inferred from the core, which of course will 
degrade the sensitivity of the experiment. If the errors are as well-behaved as those of SDSS, 
this problem can easily be resolved by the standard device of 3-a clipping. If not, then 
more complex strategies will be required. These are likely to be among the biggest practical 
problems facing the analysis, but in the absence of real data, they cannot be further analyzed 
here. 

This work was supported in part by grant AST 02-01266 from the NSF. 
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Fig. I. — The relative number of potential transiting systems F(M V ) probed for fixed plan- 
etary radius r and semi-major axis a as a function of My. The bold line applies to a uniform 
distribution of stars - to model the immediate solar neighborhood. The thin line applies to 
a thin disk - to model a search of a large portion of the Galactic disk. The dashed lines 
indicate where the two different methods for calculating the spatial density (as described in 
§ 2) overlap in each case. The distributions are arbitrarily scaled such that F(M V = 5) = 1. 
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Fig. 2. — The main plot shows the highest A% 2 values for each of 25 runs, plotted by rank, 
for each type of observing schedule. For four of the curves, a period search is conducted 
only for periods of 3.0 days to 3.1 days. Shown are a random schedule {solid), a regular 
schedule {short- dashed), a regular schedule with 14.4-minute pair-spacing {long- dashed), a 
regular schedule with 43.2-minute pair-spacing {dot-dashed). The dotted line is the random 
schedule with a period range of 3.0 to 3.5 days. The inset plot shows the main simulation 
of 1000 data sets. The 0.1% highest value is at 36.6. 



